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Progress has been made on the calculation of R 3 , the three-jet rate in e + e~ 
annihilation, in the k±_ (Durham) scheme. Using the coherent branching 
formalism J| ; [!| f|/, an explicit expression for R 3 is calculated. In this, lead- 
ing and next-to-leading large logarithms (LL and NLL) are resummed to all 
orders in QCD perturbation theory. In addition to exponentials an error 
function is involved. 
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The process e + e~ annihilation to hadronic jets has proved itself amply 
as a powerful tool for the analysis of QCD. Yet, as it is always jets which 
are observed rather than the QCD objects in which one is interested, it is 
necessary to have some mechanism for deciphering useful information from 
the observed hadronic final states. The k± or Durham algorithm is just 
one of several jet clustering algorithms which are used to define jets. As the 
name suggests it is based on the transverse momenta of the observed hadrons 
and works as follows: 

1. Define a resolution parameter y cut 

2. For every pair of hadrons hk, hi compute 

2(1 - cosd ki) /jp2 ^ 

Q 2 



VM = V ^ kl, min{ElE 2 



3. If yij is the smallest value of y^ computed in 2. and y^i < y cu t then com- 
bine the momenta Pi,Pj into a single jet ('pseudoparticle') Pij according 
to some recombination prescription. 

4. Repeat this procedure from step 2. until all pairs of particles and/or 
pseudoparticles have yui > y C ut- Whatever objects now remain are 
called jets. 

This fulfills the conditions that the algorithm should be infrared and 
collinear safe and that it should be subject to small hadronization corrections. 
Work by Catani, Webber et al.|l[ has given implicit integral expressions for 
the jet rates R n {y C ut) which will exactly sum leading and next-to-leading 
powers of \n(l/y cut ) to all orders of perturbative QCD. However, what is 
required for phenomenological applications is a set of explicit expressions. 
Numerical evaluation of the integral equations mentioned above has provided 
some useful information and has also acted as a valuable check for our results. 

The broad aim of recent work has been to find a general expression for the 
n-jet rate R n {y cu t) and this is most easily expressed in terms of the generating 
function (f>(Q, Qo; u) where 

Ql = VcutQ 2 (i) 

and u is a jet label. The n-jet rate is then 



1 ( d \ n 

Rnhjcut) = ( t^- j <f>(Q, Qo; u) 



(2) 
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From the application of the coherent branching formalism P, [|, |J, the gen- 
erating function obeys the following implicit coupled equations M: 



, Qo] u) = w 2 exp 1 2 J dqT q (Q, q) [<j) g (q, Q ; u) - 1] 



(3) 



and 



(Q,Q ;u) = uexp^y dq{V g (Q,q)[(j) g (q,Q Q ;u) - 1] 

1+u [ Q dqT f (q)exp( f" dq'{[2T q (q,q') -T g (q,q% 
JQo \ JQo 



'Qo \JQo 

{Mq',Qo;u)-i] + r f (q')} 

Where the emission probabilities are defined as 

r,(Q,<z) 
r 9 (Q,g) 



2C F a s (g) 






2C A a s (g) 






N f a s (q) 




3n q 





(4) 

(5) 
(6) 
(7) 



For R2{y C ut) one has a simple exponentiation and an explicit expression (8) for 
the resummed leading and next-to-leading logarithms [|TJ, Q is straightforward 
to obtain: 

— r — [o — L) — Po 



R 2 = exp 



12 



(8) 



where L = \n(l/y cut ), a = a s (Q)/n. Applying the formalism to the case of 
Rz{Vcut) one finds that 



R 3 = 2R 2 f Q dqT q (Q,q)exp \- f dq'(T g (q, q') + T f (q')) 

JQo I JQo 



(9) 



Substituting the emission probabilities into equation (9) - and bearing in 
mind the running of a s - gives a virtually intractable nested integral. There 
are also no obvious means by which to isolate before integrating terms which 
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will contribute to LL or NLL parts of R3. Our approach centred on perform- 
ing the calculation in two easily tractable stages. 

Firstly the integral equation was evaluated with (3q set to zero to give the 
terms in R 3 which are independent of /3q. 



AC F aR 2 
C FdL 



Qdq ( Q 3\ 

- 1 ln 7 - 4 1 ex P 



Qo q 
exp 



l' q ,1 Q 
-2C A a / dq'-\n^- 

JQo q' q' 



(10) 



CpaL 



- U,/^erf(v / I)(3(l - C F aL 2 ) - 2L) 

I Zs V i\ 



+ ^(l-expM])(A + 3 C F )} 

where A = Ca(iL 2 /4 and erf(x) is defined by 

erf(x) = — = / e~ t2 dt 

V 7T Jo 



(11) 



(12) 



Note that we have used @ = (11CU — 2iV/)/3 to write Tf in terms of (3$. 
Secondly 



dR, 



d 
Wo 

exp 



A>=0 

2R 2 



rQ 

/ dq 

JQo 



I* dq 

JQo 



2C F a a (q) LQ _ 3\ 
ixq \ q A) 
,2C A a s (qJ_^q_ _ (3 Q a s (q') 



irq' 



2-nq' 



Po=o 



with 



<*s(Q) = 



a s {jj) 



4tt pb 



(13) 



(14) 



was evaluated, thus giving the coefficients of the terms in R 3 which are pro- 
portional to (3 . This approach hinges on the assumption that R 3 at NLL 
level is linear in (3 . To justify this 



Ri 



d 2 R 3 



(15) 



A>=o 
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was evaluated and was found to contain only subleading terms (that is terms 
~ a n L 2n ~ p with p > 1), hence validating our method. 

Any subleading terms which appeared as a result of this procedure were 
discarded since they would have necessarily been incomplete. Combining the 
results of the two stages according to 



R 3 — i?o + [3 R 3 



(16) 



gives: 



CpaL 



-cxp 



C F aL 2 



-erf(VA) 



3-2L + 



(3 aL 2 ( CfclL 2 



5 H^'^ 



+^(exp[-A]-l 



C F (3 aL 2 



o 



2C> 



(17) 



This has been expressed in terms of error functions and exponentials for 
purposes of familiarity; but it should be noted that an equally succinct rep- 
resentation in terms of degenerate hypergeometric functions is possible and 
work on LL expressions for higher jet rates has shown that these functions 
appear to occur naturally. 

One might speculate on the possibility of making a change of renormal- 
isation scale such that one could absorb any terms proportional to j3o- This 
can be done for R 2 where the (3q term in equation (8) can be removed by re- 
placing a s (Q) by a s {y],!^Q). A closer look at equation (11), however, reveals 
that p enters R 3 through a combination of the emission probabilities as well 
as through the renormalisation scale; so it seems likely that no scale can be 
chosen so as to remove all (3 dependence. 

In conclusion, we have found an explicit form for R 3 {i)cut) which resums 
leading and next-to-leading large logarithms, ln(l/y cut ), to all orders in per- 
turbative QCD. Phenomenological applications of equation (17) will be dis- 
cussed in forthcoming work. 
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